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Gravitational lensing provides a strict test of cosmogonic models. 
Detailed numerical propagation of light rays through a universe hav- 
ing a distribution of inhomogeneities derived from the standard CDM 
(cold dark matter) scenario with the aid of massive, fully nonlinear 
computer simulations is used to test the model. It predicts that we 

C/} should have seen far more widely split quasar images than have been 

found. These and other inconsistencies rule out the COBE normal- 
ized CDM model with = 1 and Ho = 50km/s/Mpc; but variants of 

^ this model might be constructed, which could pass the stringent tests 

■^j- provided by strong gravitational lensing. 

Gravitational lensing directly measures fluctuations in the gravitational po- 
O tential along lines of sight to distant objects. In contrast, the conventional tools 

^ for comparing theories with observations rely on either galaxy density or veloc- 

ity information, both of which unavoidably suffer from the uncertainties with 
Q regard to density or velocity bias of galaxies over the underlying mass distribu- 

q tion, hampering our attempts to understand the more "fundamental" questions 

^} concerning the mass evolution and distribution. Thus, gravitational lensing pro- 

vides a powerful independent test of cosmogonic models 1 ' 2 . Each model for 
the development of cosmogonic structure (e.g. the HDM [hot dark matter] or 
CDM scenario) has at least one free parameter, the amplitude of the density (or 
potential) power spectrum. But now in the light of COBE observations 3 , that 
parameter is fixed by the (±15%) determination on the 5° — 10° scale in the 
linear regime. With its amplitude fixed, a secure determination of the potential 
fluctuation on any scale provides a test; any single conflict between the theory 
and reality can falsify the former. The most leverage is obtained for tests made 
on scales as far as possible from the COBE measurements. The reason is that 
all models have an assumed power spectrum that passes through the COBE nor- 
malization point at the very large comoving scales (A % lOOOMpc) fixed by that 
measurement. Since the slope of the power spectrum is a primary model depen- 
dent feature, the maximum variations amongst models occurs typically at the 
smallest scales. Thus one looks for tests at scales as small as possible, but they 
should not be so small as to be greatly influenced by the difficulty in modelling 
the physics of the gaseous, baryonic components (< lOkpc). Thus critical tests 
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are best made on scales O.OlMpc < r < IMpc. The purpose of this paper is 
to use gravitational lensing from matter distributions on these scales to test the 
standard CDM scenario. 

The model simulated here is the "standard" cold dark matter (CDM) sce- 
nario with O = 1, A = and Ho = 100/i = 50km/s/Mpc. Normalization, taken 
from the COBE first year results 3 , corresponds to erg = 1.05 4 . In order to allow 
for the existence of very large-scale waves, we first ran an L = 400/i _1 Mpc size 
box with 500 3 = 10 8 ' 1 cells and 250 3 = 10 7 ' 2 particles. In addition, in order to 
have detailed small scale information we reran a total of 10 independent simula- 
tions with L = 5/i _1 Mpc, having 500 3 = 10 8 - 1 cells and 250 3 = 10 7 - 2 particles. 
Knowing the distribution of overdensities on the 5/i _1 Mpc scale from the large 
simulation, we can statistically convolve the small and large scale runs to pro- 
duce simulated sheets or screens of matter spaced 5/i _1 Mpc apart between the 
observer at z = and a putative galaxy or quasar in the source plane at z = z s . 
A large number of independent runs (ten runs were simulated) is required so 
that identical structures do not repeat along a line of sight. The details of 
the convolution method and tests of it using a high resolution P 3 M simulation 
provided by Bertschinger and Gelb 5 will be presented in a subsequent paper. 
The method is statistically reliable for describing structures in the range 30kpc 
< ALh < 1.2Mpc, which corresponds roughly to splitting angles 5 < < 200 . 
On these scales we expect that dark matter dominates over baryons so that a 
dark matter only simulation is approximately valid. 

A very preliminary attack on this problem has been presented in Paper I 2 . 
In that work no ray tracing was done. Rather we simply checked whether or not 
mass accumulations were greater than the critical level 6 at which multiple imag- 
ing will occur. In addition to the much better method used here, the convolution 
algorithm has been modified and improved over the one adopted in Paper I. 

In our ray tracing routine we use the multiplane lens equations 7 , and speed 
the calculation of deflection angles by use of the hierarchical tree code 8 , with 
typically of order of 200-300 (grouped) screens used for each ray tracing. Ampli- 
tude on the image plane is simply given by differential area within a bundle of 
rays as compared to what it would have been had the propagation been through 
a universe with smoothly distributed matter. Figure (1) shows the amplitude 
distribution for sources of a redshift z s = 3. When a given region on the source 
plane is reached by rays in separate disconnected regions of the image plane, the 
observer would see multiple images of the same object. 

We have computed the distribution of magnifications for single and multiply- 
imaged point sources as a function of z S} multiplicity of images, and distribution 
of angular splittings. In addition, for extended sources (mock galaxies), we have 
computed the expected shape distortions, frequency and properties of the giant 
arcs that would be seen, when the source are lensed by intervening clusters 9 . Fig- 
ure 2 shows the probability of a splitting with separation of images greater than 
5 and magnitude difference less than 1.5 mag as a function of source redshift. In 
fact "amplification bias" 6 , 10 ' 11 will increase the probabilities over those shown 
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Probability of a splitting with separation of images greater than 5 and magnitude difference 
less than 1.5 as a function of source redshift. 



in Figure 2 by a significant amount. Splittings larger than 5 should be common 
(when several thousand quasars have been examined), if this cosmogonic model 
were correct. 

Probably the single most revealing statistic is the distribution of image sep- 
arations expected for multiple sources as shown in Figure 3 for z s = 1, 2 and 
3. Notice that very large splittings should be the rule. Also revealing is the 
distribution of expected lens redshifts as shown in Figure 4. The lenses them- 
selves should be close enough to be seen in almost all cases. On this issue the 
recent observation 12 of a lens candidate for the double quasar QSO2345+007 is 
extremely relevant. The separation of the two images is 7 .06, the quasar red- 
shift is z s = 2.15 and the putative lens is at z\ = 1.49. We see from Figure 
4 that, although 7 arcsecs separation can be produced in the CDM model, the 
probability that the lens is as far away as z = 1.49 is very small (2%) due to the 
relatively late formation of structure in this model. In open models structure 
formation occurs earlier. 

It appears that all three of these results (shown in Figures 2,3,4) are seri- 
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Figure 3 

Multiple-lensing probability distribution as a function of image separations for sources at Z s = 
1, 2 and 3. Also shown as long dashed curve is the observed distribution 13 . 

ously in conflict with the existing observations. In particular, we find that the 
standard CDM model predicts that 0.0007 of all lines of sight to z s = 1, 0.0014 of 
all lines of sight to z s = 2 and 0.0020 of those to z s = 3 will be multiply-imaged 
with angular splittings >10" and amplification ratios of less than 1.5 magnitudes. 
Surveys 12-18 and occasional serendipitous discoveries have revealed 27 confirmed 
or possible multiply-imaged QSO's according to a recent compilation 19 . Detailed 
analysis 20 ' 21 of the most statistically useful of these surveys 14 ' 15 ' 16 yields a lens- 
ing rate in the vicinity of a few tenth to one percent, consistent with the CDM 
predictions quoted above making allowance for plausible magnification biases 22 . 
However, as shown in Figure 3, all observed QSO lens systems have image split- 
tings of less than 10 , and the large majority, less than 5 . 

This sharply contradicts and thus falsifies the model. Since the large split- 
ting, modest brightness ratio systems predicted by the model would be typically 
much easier to detect and recognize than those 27 which have actually been 
found, no escape by appeal to observational selection seems possible. 

This failing of the model is not presented as an entirely new result, but 
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Integrated lensing probability distribution as a function of expected lens redshift. The X symbol 



indicates the recent observation of of a lens candidate for the double quasar QSO2345+007. 



only as a new and more robust manifestation of a previously recognized prob- 
lem, namely the excessively deep potential wells produced by the dark matter 
component in COBE normalized standard CDM 23 ' 24 ' 25 . These excessively deep 
potential wells lead both to excess galaxy pair-wise velocity dispersions and to 
the predicted excessive rate of large splitting lensing events. The virtue of the 
lensing test is that it is independent of other tests and is not subject to the same 
caveats concerning "bias" of galaxies with respect to dark matter. 

Are there variant models that would not fail these tests? The simplest 
change that one can make in the scenario is to reduce the density of the clumped 
material in the universe, as it is the large number of large mass concentrations 
that produce the over- abundance of large splittings. Thus a lower is clearly 
useful, but it would be premature to argue that the present results by themselves 
indicate < 1, as many other properties of the scenario ("temperature" of the 
dark matter, shape of the power spectrum etc) contribute to lensing properties. 
However the directness of gravitational lensing as a test for the growth of inhomo- 
geneities, coupled with the rapidly increasing power of computers and numerical 
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algorithms, makes one optimistic that calculations of the type reported on here 
should become a major tool for testing and discriminating among competing 
cosmological scenarios. 
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